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PREFACE 


The  Rand  Corporation,  under  sponsorship  of  the  Tactical  Technology 
Office  of  the  Defense  Advanced  Research  Projects  Agency,  has  been  en- 
gaged in  analysis,  development,  and  evaluation  of  hydrodynamic  design 
criteria  for  submersible  vehicles.  Boundary-layer  control  including 
shaping,  pressure  gradients,  suction,  and  heating  has  been  the  main 
technological  basis  for  these  investigations. 

Successful  boundary-layer  control  depends  on  appropriate  manipula- 
tion of  the  laminar  boundary  layer  to  maintain  stability  and  delay  tran- 
sition to  turbulence.  Thus,  the  effects  of  surface  heating  and  pressure 
gradients  on  properties  of  the  laminar  boundary  layer  in  water  are  cen- 
tral to  progress  in  low-drag  technology. 

This  report  provides  a simplified  and  accurate  approximate  method 
I',  calculating  the  properties  of  heated  laminar  boundary  layers  in 
water.  When  coupled  with  appropriate  stability  or  transition  criteria, 
it  provides  a systematic  and  efficient  approach  to  the  optimization  of 
low-drag  design. 

The  method  presented  here  should  simplify  study  of  the  interplay 
between  pressure  gradients  and  surface  heating  on  laminar  characteris- 
_ics  to  the  point  where  a hand  calculator  can  be  used  to  compute  laminar- 
flow  characteristics  for  the  case  of  cjnstant  wall  temperature.  Addi- 
tional work,  still  in  progress,  will  extend  this  approach  to  the  case 
of  variable  wall  temperature. 

This  report  should  be  useful  to  hydrodynamicists , to  designers  of 
submersibles,  and  to  others  interested  in  applyinp  fluid  mechanics  to 
improve  the  performance  of  underwater  vehicles. 

Other  related  Rand  publications  are; 

R-1752-ARPA/ONR,  Low-Speed  Boundary -Layer  Traveition  Workshop^ 

W.  S.  King,  June  1975. 

r_1789_ARPa,  Controlling  the  Separation  of  Laminar  Boundary 

Layers  in  Mter:  Heating  and  Suction^  J.  Aroesty  and  S.  A. 
Berger,  September  1975. 
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R-1863-ARPA,  The  Effects  of  Wall  Tempevature  and  Suction  on 
haminar  Boimdxry-Layer  Stability,  W.  S.  King,  April  1976. 

9 

R-189<'-ARPA,  "e  Stability  Theory  and  Boundary-Layer  Transi- 
tion, S.  A.  Berger  and  J.  Aroeaty,  February  1977. 

R-1907-ARPA,  Buoyancy  Cross-Flow  Effects  on  the  Boundary  Layer 
of  a Heated  Horizontal  Cylinder,  L.  S.  Yao  and  I.  Catton, 
April  1976. 

R-1966-ARPA,  The  Buoyancy  and  Variable  Viscosity  Effects  on  a 
Water  Laminar  Boundary  Layer  Along  a Heated  Longitudinal 
Horizontal  Cylinder,  L.  S.  Yao  and  I.  Catton,  February 
1977. 

R-2111-ARPA,  Entry  Flow  in  a Heated  Tube,  L.  S.  Yao,  June  1977. 

R-2164-ARPA,  The  Effects  of  Unsteady  Potential  Flow  on  Heated 
Laminar  Boundary  Layers  in  Water:  Flow  Properties  and 
Stability,  W.  S.  King,  J.  Aroesty,  L.  S.  Yao,  and  W. 
Matyakiela,  in  process. 


The  authors  wish  to  express  their  appreciation  to  C.  Gazley,  Jr. 
for  reading  the  manuscript  and  offering  several  helpful  suggestions. 

The  authors  also  wish  to  thank  W.  H.  Krase  and  L.  S.  Yao  for  their  care- 


ful reviews. 
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Thwaitos^’^  and  otliers^^^  have  used  integral  methods  to  compute 
parameters  of  laminar  boundary  lavers  with  constant  fluid  properties. 

These  methods  use  simple  correlations  of  universal  parameters  based 
principally  on  similar  boundary  layers.  The  displacement  thickness, 
momentum  thickness,  and  wall  shear  stress  can  be  computed  by  these 
metlx'ds  for  nonsimilar  boundary  layers. 

This  report  is  concerned  with  the  extension  and  application  of 
these  methods  to  heated  water  boundary  layers.  The  viscosity  of  water 
is  strongly  temperature  dependent,  and  thus  variable  fluid  property 
effects  are  significant  for  heated  water  boundary  layers.  The  Thwaites 
integral  method  has  been  extended  to  include  variable  fluid  properties. 
Numerical  solutions  of  heated  water  wedge  flows  arc  tabulated  for  a 
range  of  wedge  angles  and  temperature  differences.  Universal  param- 
eters are  correlated  from  these  solutions  for  use  with  the  integral 
method.  The  method  has  been  tested  for  the  nonsimilar  Howarth  retarded 
flow  with  t^  = 10A°F  and  t^  = 32°F.  Errors  of  only  1 or  2 percent  are 
found  for  displacement  thickness,  momentum  thickness,  and  wall  shear 
stress.  This  is  a particularly  stringent  test  because  of  the  adverse 
pressure  gradient  in  the  Howarth  flow,  where  heating  causes  the  separ- 
ation point  to  move  30  percent  further  down  the  surface. 

The  curvature  of  the  velocity  profile  at  the  surface,  f"'(0),  de- 
pends on  the  viscosity  gradient  at  the  surface,  which  In  turn  is  deter- 
mined by  the  heat  flux  at  the  surface.  These  three  additional  parameters 
can  be  computed  once  the  Nusselt  number  is  known.  The  Lighthill  high- 
Prandtl -number  approximation  has  been  extended  to  permit  computation 
of  the  Nusselt  number  for  boundary  layers  wltli  variable  fluid  proper- 
ties. Nusselt  numbers  computed  with  this  approximate  method  for  the 
heated  Howarth  flow  are  reasonably  close  to  the  exact  numerical  solu- 
tions for  that  flow,  except  in  the  region  near  separation. 

Heated  water  boundary  layers  have  been  further  explored  with  simpli- 
fied shear  modeling,  such  as  Couette  flow.  The  Couette  flow  model  works 
very  well  at  stagnation  points,  where  the  boundary  layer  does  not  grow. 
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NOMENCLATURE 


a,b  = constants 


c = pp/Pg  Pg 

pk  1 
“ (pk)^  Pr^ 

c = heat  capacity 
P 

f = dimensionless  stream  function,  f “ u/U 
F - universal  parameter,  2T  - 2(2  + H)X 
H * universal  parameter,  <5*/9 

j = 1 for  axisymmetric , 0 for  two-dimensional  bodies 
k = thermal  conductivity 
L - an  arbitrary  length 


m - (x/U)  dU/dx 
Nu  = Nusselt  number,  hx/k^ 

Pr  * Prandtl  number 

r * radial  distance  from  the  symmetric  axis 
Re^  = Reynolds  number,  x 
t = temperature 

T = universal  parameter,  ^/pU 

u * X component  of  velocity 

U « velocity  outside  the  boundary  layer 


V 

X 

y 

3 


y component  of  velocity 
coordinate  along  the  surface 
coordinate  normal  to  the  surface 


2m/ (m  + 1) 


6 
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6*  = displacement  thickness,  / (1  - u/U)  dy 

•'o 

n = sl...llarlty  variable,  y/6(x) 

6 “ momentum  thickness,  / (1  - dy 

*/q  U U 

0 - e/6 

X - (e^/v  ) dU/dx 

w 

U ■=  viscosity 
V = kinematic  viscosity 
p - density 
T • shear  stress 

Subscripts 

e = at  ambient  conditions 
Iso  = Isothermal 

T * of  the  thermal  boundary  layer 
w = at  the  surface  conditions 
ref  = at  A91.69®R  (32”F) 


1 


INTRODUCTION 


The  key  to  boimdnry-layer  roiitrol  is  ttie  maintenance  of  laminar 
velocity  profiles  that  promote  hydrodynamic  stability  and  delay  separ- 
ation. I’or  undersea  app  1 icat ions , surface  heating  is  a possible  means 
of  achieving  this  control . ^ Heating  promote.s  stability  through 

the  interplay  among  the  thermal  boundary  layer,  the  temperature-dependent 
viscosity  of  water,  and  the  momentum  balance  in  the  crucial  region  near 
the  wall. 

To  take  full  advantage  of  these  phenomena,  the  surface  heating  dis- 
tribution should  be  optimized,  so  that  total  heat  required  is  minimized 
subject  to  realistic  constraints.  As  part  of  a study  to  simplify  opti- 
mization and  design  methods,  we  have  developed  an  approximate  method 
for  determining  laminar  boundary-layer  characteristics  for  cases  of 
constant  surface  overheat.  The  extension  of  this  approach  to  the  case 
of  variable  surface  temperature  is  currently  under  investigation  at  Rand. 

Approximate  methods  are  useful  in  optimization  studies  because  nu- 
merical integration  of  the  boundary-layer  equations  is  unnecessary.  If, 
for  example,  the  optimal  heating  distribution  for  a given  body  shape 
is  required,  it  may  be  necessary  to  perform  hundreds  of  numerical  calcu- 
lations of  the  boundary-layer  equations,  with  consequent  cost  in  computer 
time.  Approximate  methods  can  reduce  this  cost  considerably  and  can 
also  yield  better  intuitive  understanding. 

For  a method  to  be  useful  in  predicting  stability,  transition,  and 
separation,  it  must  supply  information  about  the  velocity  profiles  that 
influence  these  phenomena.  For  example,  the  critical  Reynolds  number 
for  stability  of  Tollmein-Schlichtlng  waves  has  been  shown  to  correlate 
with  the  shape  factor,  H * 6*/0,  and  to  depend  on  the  curvature  of  the 
velocity  profile  at  low  values  of  Integral  methods  can  be  used 

to  determine  these  parameters  with  sufficient  accuracy  for  use  in  opti- 
mization studies. 

We  have  extended  Thwaites's  integral  method  to  ^ae  case  of  water 
flowing  over  constant  temperature  surfaces.  Thwaites'u  method  depends 


*Although  this  method'  is  conventionally  called  after  Thwaites, 
Walz,  Thwaites,  and  others  were  responsible  for  its  development.  See 
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'in  t'xti'iiKivt’  rorrr  1 at  idiis  of  parameters  derived  from  numerically  exact 
l)('imdarv- 1 <ivc-r  solutions,  raLher  than  a particular  assumed  form  of 
vi-loi'ity  profile.  The  key  simplification  in  the  Thwaites  nuthod  arises 
from  tlu’  ohservathin  that  the  parameter  F = 2T  - 2(2  + 11)A  c.in  be 
approx  im.'it  el  V expressc'd  as  a - hX,  wlu're  a aad  b ar  constants  derived 
from  a best  fit  to  a graph  of  F vs  X from  a computed  series  of  different 
pressure  gradients.  In  the  constant  property  (ase,  a = .44  and  b = “i . 38 
.ire  widely  used  v.iliies.  We  have  found  that  an  approximate  linear  rela- 
tion between  F and  A is  still  valid,  nut  the  constants  a and  b now  de- 
]-)end  on  wall  overlie, at  aiul  ainbiriit  t (’mper.a  t ure . As  in  Tliwaites's  origin, il 
method,  the  momentum  thickness  il  i s t r i bu  t i on  can  be  obtained  from  a simple' 
(piadratui'e,  while  shear  stress,  shape  factors,  and  profile*  enrvaturc  ean 
be*  eibtained  from  ediarts  or  t.ahles. 

Tlu'  Thwaites  eipproach  is  a simple,  ee-oneimi  c.i  1 , and  practie-al  scheme 
fell*  the  c.i  1 cu  1 ,it  iem  of  nt)nsiniilar  boundarv  layers  ,ind  is  .ae-eiirate  ene>ugh 
f eir  most  applications  except  in  the  neighbeir hood  eif  laminar  separation 
With  the  extension  we  present  here,  it  is  possj  tei  compute  leiminar 
boundary  e harne' t er  i st  i r s e>f  water  flowing  over  e'onstant  temperature'  siir- 
fae-es  with  a hand  cale-ulator  .alone. 

Ref.  2 for  a historiceil  survey.  The  method  we  adopt  here,  because  of 
its  reliance  on  similarity  solutions,  is  e leisest  to  that  proposed  by 
Wa 1 X for  the  constant  property  ease. 
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II.  INTEGRAL  METHOD 


The  conservation  equations  for  mass,  momentum,  and  thermal  enoi 
in  a steady  boundary  layer  are 


(la) 


(lb) 


(Ic) 


where  j = 1 for  axisymmetric  flow  and  j = 0 for  two-dimensional  flow. 
If  the  density  is  not  a strong  function  of  temperature  (as  for  water), 
the  momentum  equation  can  be  integrated  across  the  boundary  layer, 
yielding 


^w  ^ 26  + 5*  ^ ^ ^ (2a) 

2 “ dx  U dx  r dx  ’ 

pU 


where 


e = 


(2b) 


and 


6* 


(2c) 


f 1 2^ 

The  integral  method  developed  by  Thwaites'  ’ has  been  useful  for 
unheated  boundary  layers.  In  this  method,  universal  parameters  corre- 
lated from  a number  of  known  solutions  for  similar  and  nonsimilar  flows 
are  applied  to  boundary- layer  computations  for  any  body  shape.  The  con- 
venient universal  parameters  F,  H,  and  T,  defined  as 
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F = 2T  - 2(2  + H)A 
H >=  6*/e 
T = T 0/yU  , 

W 

are  correlated  as  functions  of  X,  defined  as 

, _ ^ dU 
V dx 


(3a) 

(3b) 

(3c) 


(4) 


The  Integrated  momentum  equation  (2a)  can  be  expressed  as 


1 ^ / r^^X  \ ^ F(X) 
^2j  dx  \ dU/dx / U 


If  F Is  linear  In  X,  F ■=  a - bX  (as  shown  to  be  approximately  the  case  by 

(1  2) 

Thwaltes  ’ for  constant  fluid  properties),  then  Eq.  (5)  can  be  formally 
integrated,  yielding 


For  a given  U = U(x),  9 can  be  obtained  from  this  expression;  X is  then  ob- 
tained from  its  definition  and  H and  T are  then  found  from  their  correla- 
tions as  functions  of  X.  With  this  procedure,  6(x),  6*(x),  and  t (x)  can 

w 

be  calculated  for  flow  over  a surface  of  any  shape,  once  U(x)  is  known.  It 

(1  2) 

has  been  shown  ’ that  the  Thwaites  integral  method  accurately  (with- 
in about  5 percent)  determines  H and  skin  friction  for  general  nonslml- 
lar  laminar  boundary  layers  with  constant  fluid  properties,  except  in  the 
region  near  separation. 

The  variable  viscosity  analogs  to  Eqs.  (3)  through  (6)  are  almost 
identical  to  these  constant  fluid  property  equations.  Equations  (2)  hold 
when  the  viscosity  is  temperature  dependent.  If  the  definitions  in  Eqs.  (3) 
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and  (4)  are  modified  only  to  the  extent  of  evaluating  v and  v at  the  wall 
temperature  (replace  y and  v with  y^  and  v^) , then  Eq.  (5)  holds  for 
variable  viscosity.  The  simple  algebraic  derivation  of  Eq.  (5)  for  the 
variable  viscosity  case  is  nearly  Identical  to  that  for  the  constant 
fluid  properties  case.  Again,  if  F is  assumed  linear  in  X,  Eq.  (6)  will 
hold  (with  V replaced  by  v^).  In  fact,  for  a given  constant  wall-to- 
ambient-teraperature  difference,  it  will  shortly  be  shown  that  F is  indeed 
approximately  linear  in  X.  Thus,  Eq.  (6)  for  6 (x)  holds  for  heated  water 
boundary  layers  as  well  as  for  unheated  ones,  but  the  universal  functions 
F(X),  H(X),  and  T(X)  will  depend  on  the  wall  and  ambient  temperatures  as 
well  as  on  X.  These  correlations  can  be  obtained  from  numerical  solutions 

for  wedge  Clows. 

For  laminar  boundary  layer  flow  of  a fluid  with  temperature-dependent 
fluid  properties  over  a wedge  (j  = 0)  or  a cone  (j  = 1)  with  a constant 
surface-to-ambient-temperature  difference,  the  boundary  layer  equations 
(la  through  c)  reduce  to  two  ordinary  differential  equations,  similarity 
equations.  Thus,  these  wedge  and  cone  flows  are  convenient  for  correlat- 
ing the  universal  functions  F,  H,  and  T.  The  continuity  equation  (la) 
can  be  eliminated  by  defining  a stream  function  as 


r^pu 


3y 


r^  pv 


3x 


(7) 


Now  if  the  Mangle r-Levy-Lees  transformation  from  x and  y coordinates  to 
S and  n coordinates  is  introduced, 

dC  = p y U(r/L)^^  dx 

6 6 

dn  = [pU/(2C)*^](r/L)^  dy  , (8b) 

where  L is  any  arbitrary  length,  and  a dimensionless  stream  function,  f, 
is  defined  such  that 


<Kx,y)  - (2^)^  L^f(5,Tl)  , 

(8) 


(9) 


then  the  boundary  layer  equations  (1)  become 
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(cf")'  + ff" 


(10a) 


where 


(V')  fS'  - 2c(f  If  - g'  If  ) 


= 2i  ^ 

■”  U dC 


_ Pk  1 

‘^k  " (pk)  Pr 
e e 


(lOb) 


(lla) 

(llb) 

(llc) 

(lid) 


for  any  constant  T^,  and  where  prlr.es  denote  derivatives  with  respect  to  n- 
The  fluid  property  ratios  are  functions  of  local  temperature,  or  of  g 
only.  If  6 and  all  boundary  conditions  are  independent  of  then  Eqs.  (10) 
become  independent  of  C*  For  flow  over  either  a wedge  or  a cone,  the 
velocity  outside  the  boundary  layer  is  of  the  form^^^  U = ^ wedge 

or  a cone,  ^ is  also  proportional  to  a power  of  x.  Thus,  for  a wedge  or 
a cone,  B (from  Eq.  (11a))  is  a constant.  For  an  isothermal  impermeable 
surface  without  slip,  the  boundary  conditions  at  the  surface  (n  = 0)  are 


f(0)  = f*(0)  = 0 g(0)  = 0 


(12a) 


if  we  set  T_  - T . To  match  the  ambient  temperature  and  velocity  at  the 
0 w 

edge  of  the  boundary  layer,  the  boundary  conditions  as  n “ are 

f'(“)->l  g(“)  1 . (12b) 

Therefore,  for  an  isothermal  wedge  or  cone,  B and  the  boundary  conditions 
are  in  fact  independent  of  K,  and  so  Eqs.  (10)  reduce  to  the  similarity 
equations 
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(cf")'  + ff"  + e 


- - (f)^ 


(C^g’)'  + fg' 


= 0 
= 0 .• 


(13a) 

(13b) 


The  transformed  velocity  profiles  for  cones  are  the  same  as  those  for 

wedges;  we  need  only  correlate  F,  H,  and  T with  wedge  flows. 

Numerical  solutions  of  wedge  flows  (Eqs.  (13))  with  water  have  been 

obtained  for  various  values  of  wedge  angle  and  wall-to-ambient-temperature 

difference  for  an  ambient  temperature  of  67°F  (see  Table  1).  The  water 

(9) 

property  correlations  used  were  tliose  of  Kaups  and  Smith  (see  the  appen- 
dix), but  density  was  taken  as  constant. 

The  above  transformed  boundary  layer  equations  are  convenient  for 
numerical  study,  but  the  true  magnitude  of  the  heating  and  pressure 
gradient  effects  on  the  shear  stress,  for  example,  may  not  be  obvious 
from  the  transformed  shear  stress,  f"(0).  To  explore  the  combined  ef- 
fects of  pressure  gradient  and  heating,  and  to  offer  further  physical 
insight  into  these  effects,  the  following  simplified  shear  models  have 
been  developed.  The  simplest  is  to  assume  that  heating  does  not 
change  the  shear  in  the  boundary  layer,  the  "isothermal  shear  model. 

Since  r is  proportional  to  c f"(0), 
w ^ 


t"(0)  „ ,-l  (14) 

” " ■ 

Another  shear  model,  the  Couette  flow  model,  assumes  that  the  shear 
stress  (but  not  the  velocity  gradient)  is  a constant  for  the  entire  vis- 
cous boundary  layer.  The  Couette  flow  model  requires  that  some  distribu 
tion  of  viscosity  be  assumed.  We  have  assumed  that  l/p  is  distributed 
linearly  from  l/p^  at  the  wall  to  l/p^  at  the  thermal  boundary  layer 
edge.  This  viscosity  distribution  assumption  corresponds  closely  to  a 
linear  temperature  distribution.  The  Couette  flow  model  also  requires 
the  assumption  that  the  viscous  boundary  layer  thickness  is  not  changed 
appreciably  by  heating;  this  can  be  shown  to  be  the  case  from  the  top 
part  of  Table  1.  This  model  results  in 


Table  1 


WEDGE  FLOW  NUMERICAL  SOLUTIONS  WITH  WATER  AT 

AN  AMITIEN!  TEMPERATURE  t = b7°F 

o 


f"(0),  0,  and  H as  functions  of  6 and  At 


A - 20  F 
c . 0.7758 


- 30  K 
r - 0.6927 


f'  (0)  and  Nii/Re  as  functions  of  B and  At 


-0.18  0.1907 

-0.15  0.1459 

-0.10  0.07470 

-0.05  0.005474 

-0.025  -0.02865 
0.0  -0.06256 

0.025  -0.09626 
0.05  !-0. 1298 
0.10  -0.1963 
0.20  -0.3280 
0.50  -0.7152 
1.00  -1.147 


Nu/Re 

1 „ . 

f"(Q) 

U 

Nu/Re 

0.4860 

0.1983 

0.5017 

0. 5359 

0. 1382 

0.5489 

0.5894 

0.04500 

0.6010 

0.6291 

-0.04420 

0.6401 

0.6463 

-0.0H782 

0.65' 

0.6623 

-0.1310 

0.673 

0.6774 

-0. 1737 

0.6881 

0. 6918 

-0.2161 

0.7023 

0.7189 

-0.2996 

0.7293 

0.7690 

-0.4640 

0.7794 

0.9110 

-0.9417 

0.9214 

-1.709 

1.212 

0.1272 

0.01141 


0.5163 
j 0.5613 
0.6121 

0. 6507 
■1.6676 

1.  6833 
0.6982 
0. 7124 
0.7393 
0.7892 

j 0.9313 


0.2023 

0.09638 

-0.06525 

-0.2169 

-0.2902 

-0.3622 

■0.4329 

■0.5029 

■0.6401 

■0.9067 

•1.665 

■2.854 
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, (15) 

1 + ^(6  /6)(c"^  - 1) 

Z i w 


where  the  ratio  of  boundary-layer  thicknesses  can  be  approximated  by 


_-l/3 

__L  ~ pr 
& 


(16) 


These  simplified  shear  models  suggest  that  heating  may  affect  f"(0) 

by  approximately  the  same  factor  for  all  pressure  gradients;  this  f"(0) 

ratio  is  shown  in  Table  2 for  wedge  flows  with  various  values  of  B,  and 

in  fact,  this  ratio  is  fairly  constant  for  a given  At.  Note  that 

the  Couette  flow  model  is  very  good  at  6 = 1,  stagnation  point  flow, 

where  the  boundary  layer  does  not  grow  (see  Table  2).  The  Isothermal 

shear  model  corresponds  closely  to  B = -0.1  (see  Table  2),  for  which  it 

can  be  seen  from  Table  1 that  the  shear  stress  (proportional  to  c^f"(0)) 

is  nearly  the  same  for  all  heating  levels. 

Returning  now  to  our  original  analysis,  we  note  that  Eqs.  (3)  for 

F,  H,  and  T can  be  expressed  in  terms  of  the  transformed  variables. 

With  the  density  taken  as  constant,  for  a wedge  transformation  (8b)  be- 

u 

comes  n = y/6(x),  where  6(x)  = /T - 8 Re""^.  Equations  (3)  become 

F = 2(1  - B)9^c"^  (17a) 

w 

H = 6*/e  (17b) 

T - 0f"(O)  (17c) 


for  isothermal  wedge  flows,  where  9 * 9/(5.  For  wedge  flows,  X can  be 
found  from  the  following  boundary  condition  at  the  wall 


(18) 


from  which 
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Table  2 

COMPARISON  OF  SIMPLIFIED  SHEAR  MODELS  WITH  NUMERICAL  SOLUTIONS 


OF 

f"(0)  FOR 

HEATED  WATER  BOUNDARY 

LAYERS  OVER  WEDGES 

(t 

= 67°V) 
e 

f"(0)/f" 

Isi 

(0) 

B 

fV  (0) 

ISO 

At  = 10°F 

1 

At  = 20°F 

At  = 30”f 

At  = 50°F 

- . 18 

. 1286 

1 . 246 

1.506 

1.778 

2.349 

- . 15 

.2164 

1.169 

1.347 

! 1.530 

1.909 

i 

- .10 

1 .3193 

1.140 

1.284 

1.432 

1.731 

- .05 

. 4003 

1.129 

1.260 

! 1 . 394 

1.662 

0.00 

.4696 

1 .123 

1.247 

1 

1.372 

1.623 

. 10 

. 5870 

1.116 

1.232 

1.348 

1.579 

. 20 

.6867 

I.  112 

1.223  ‘ 

1.334 

1.553 

.50 

.9277 

1.104 

1.208 

1.311 

1.511 

1.00 

1.233 

j 1.098 

* 1.196 

1.292 

1.477 

Isot-hermal  shear 

model 

1.141 

[ 

1.289 

1.444 

1.771 

Couette 

flow 

model 

1.100 

1.195 

1.286 

1.453 

or 


(19a) 


(19b) 


(19c) 


F,  H,  and  T have  beer  calculated  from  Eqs.  (17)  using  Table  1,  and  are 

plotted  in  Figs.  1,  2,  and  3 as  functions  of  X,  computed  from  Eq.  (19c). 

Actually,  as  indicated,  the  F and  H curves  are  plotted  as  functions  of 

X"  = c X instead  of  X,  and  F = c F is  plotted  Instead  of  F.  This  removes 
w w 

c from  F and  X (see  Eqs.  (17a)  and  (19c))  and  makes  the  resulting  curves 
w 

nearly  parallel  and  closer  together.  Fortunately,  as  is  seen  in  Fig.  1, 

F is  nearly  linear  in  X for  a given  At,  so  that  the  extension  of  Thwaites's 
method  is  straightforward.  The  form  F = a - bX  can  be  used,  but  the  con- 
stants a and  b now  depend  on  temperature.  As  shown  in  Figs.  1,  2,  and  3, 
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F = a - 

a = [0.441  - 0.0953(c"^  - l)]c"^ 

w w 

b = 5.38  + 0.415(c"^  - 1) 
w 

H = }K(A)  - 0.410(c~^  - 1) 
u w 

T = T_(^)  + (0.149  - 0.730A)(c~^  - 1)  , 

U W 


(20a) 


(20b) 

(20c) 


where  and  are  the  isothermal  H and  T values.  These  correlations 

are  based  on  wedge  flow  solutions  with  t^  = 67°F,  and  they  have  been  checked 

at  several  other  ambient  temperatures.  Perhaps  an  additional  correlating 

parameter,  such  as  Pr  or  c',  is  needed  in  addition  to  c"^.  At  t = 47°F 

e w we 

and  At  = 20'’F,  a,  b,  and  H are  predicted  by  the  above  correlations  within 
1 percent;  for  this  case,  T is  predicted  within  4 percent  except  for 
values  of  X near  separation  where  T can  be  in  error  by  as  much  as  12  per- 
cent. Hov;ever,  at  t^  = 32°F,  the  above  correlations  are  poor,  as  might 
be  expected. 
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III.  HEAT  IRANSFjfi 


Water  has  a high  Prandtl  number:  in  the  temperature  range  of  inter- 
est for  underwater  vehicles,  the  Prandtl  number  of  water  ranges  from  A^to 
13.  The  Lighthill  high-Prandtl-number  approximation  for  heat  transfer 
works  well  In  this  range;  in  fact,  for  a flat  plate  with  zero  pressure 
gradient  at  a Prandtl  number  of  one,  the  high-Prandtl-number  approxima- 
tion is  in  error  by  only  2 percent.  For  constant  fluid  properties,  the 

(10,71 

high-Prandtl-number  approximation  for  heat  transfer  is 


Nu  = 0.538A  Pr 

X 


1/3  X 


K \ r 

» I lo  'J  '' 


dx  I 

^ I 


-1/3 


(21) 


For  this  same  formula  to  hold  for  water,  with  its  temperature-dependent 
viscosity,  a reference  temperature  must  be  used.  As  demonstrated  in  the 
following  section,  this  reference  temperature  approach  meets  with  accept- 
able but  limited  success. 

The  Nusselt  number  is  a dimensionless  surface  heat  flux  or  a dimen- 
sionless temperature  gradient  at  the  surface. 


where 


Nu  = 


q X 
w 


q = k 

•».t  T. 


: (t 

- t ) 

w w 

e 

It 

In 

"w  3n 

3y 

w 

w 


When  the  Nusselt  number  has  been  determined,  the  viscosity  ratio  gradient 
at  the  wall  can  be  computed  from  the  expression 


w 


V. 


dp 

dn 


dp 

dt 


t=t 


dt 

dn 


w 


n*0 


(22) 


It  is  sometimes  useful  or  necessary  to  have  expressions  for  higher  deriva^ 
tives  of  the  velocity  at  the  wall  (e.g.,  see  comments  in  Sec.  I).  By 
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taking  the  l imit  of  Eq.  (10a)  as  n "*■  0,  f "'(0)  c ii  be  expressed  as 

f"'(0)  = -[6  +c;f"(0)]/c  , . (23) 

W w 

Higher  derivatives  of  f and  t at  the  wall  can  he  similarly  found  such 
that,  e.g.,  series  expansions  for  the  temperature  and  velocity  profiles 
near  the  wall  can  be  obtained. 

Tbi’s  0,  6*,  f " (0) , Nu  , c',  f"'(0),  and  other  parameters  can  all  be 

X w 

found  for  nonslmllar  heated  water  boundary  layers  by  simple  integrals  and 
a table  or  chart.  The  procedure,  for  an  arbitrary  pressure  gradient  of 
the  nonslmllar  type,  is  to  integrate  Eq.  (6)  to  obtain  the  momentum  thick- 
ness 0,  then  use  Eq.  (4)  and  Figs.  1,  2,  and  3 to  obtain  the  shape  factor 

H,  and  the  shear  parameter  T.  With  T and  0 known,  then  x (x)  = p UT(x)/0(x) 

w w 

is  computed,  and  the  integration  in  the  high-Prandtl -number  heat-transfer 
formula  (Eq.  (21))  csn  be  performed. 
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IV.  EXAlfl’LE:  HOWARTH'S  LINEARLY  RETARDED  H.OW 


The  methods  descrihed  In  the  previous  sections  have  been  applied 

to  a demanding  test  case  with  an  adverse  pressure  gradient  and  strongly 

(2)  , 

varying  fluid  properties.  An  exact  numerical  solution  of  Hov;arth's 

retarded  flow  [U  = u^(]  - x/8L)]  for  water  with  t^  = 104°F  and  t^  = 32°F 

has  been  compared  with  the  integral  method  solution  to  the  same  problem. 

The  correlations  (Eqs.  (20))  cannot  be  extrapolated  to  this  low  an 

amhien’:  temperature,  so  numerical  solutions  for  water  wedge  flows  were 

made  with  t = 104°F  and  t = 32“F.  F,  H,  and  T are  plotted  as  func- 
_ w 0 

tions  of  \ in  Fig.  4.  Again,  F(X)  is  a nearly  straight  line  with  a = 0.904 
and  b = 6.05.  The  comparison  between  the  exact  and  integral  solutions  (see 
Table  3)  shows  excellent  agreement  for  6 and  good  agreement  for  6*  and 
f"(0).  Note  that  8 = -.18  corresponds  to  X s .13  for  wedge  flows,  and  our 
computations  beyond  this  point  (x/L  greater  than  0.8  for  Howarth  s flow) 
are  extrapolations.  Howarth' s flow,  having  an  adverse  pressure  gradient, 
is  strongly  affected  by  variable  fluid  properties.  With  the  above  tempera- 
ture boundary  conditions,  separation  occurs  at  x/L  = 1.246,  while  in  the 
constant  fluid  properties  case  (t  = t ) separation  occurs  at  x/L  = 0.958 
(the  "8"  in  the  U equation  was  arbitrarily  chosen  to  make  x/L  near  one  at 
separation).  This  shift  of  the  separation  point  is  a dramatic  variable- 
fluid  property  effect.  Thus,  this  flow  is  a good  check  case  for  the  inte- 
gral method. 

The  hlgh-Prandtl  number  approximation  was  applied  to  this  flow,  with 
the  same  surface  and  ambient  temperatures,  using  the  numerical  solutions 
for  f"(0)  (see  Table  3),  the  average  Prandtl  number,  and  the  surface  values 
of  T and  V.  The  Nusselt  numbers  from  this  approximate  formula  are  compared 
with  exact  numerical  solutions in  Table  4.  The  difference  is  only  about 
4 percent  for  most  of  the  plate,  but  the  difference  is  much  larger  near 
separation . 


lab  If  3 


COMPAKISl  N OV  THF.  INTbCRAL  MF.TllOi)  TO  NlIMKIUCAb  SOLUT 


Hr.TARDFI)  FLOV^  (WATFR  WITH  t = 1()^°F  AND  I 


IONS  I'OK  lIOWARTH's 
= 32 '’F) 


x/1. 

Numerical  Solution^ 

Integral  Method 

P 

6Re'^/x 

X 

•k 

h /O 

f"  (0) 

ORe^/x 

X 

k 

6 /O 

f"(0) 

0.0 

0. 

0. 

.573 

2.13 

1.117 

.574 

2.12 

1.113 

. 2 

-.0247 

0.0526 

592 

2.  17 

1.022 

.593 

/.17 

1.022 

.u 

-.0544 

-.111 

.613 

2.22 

.914 

.614  1 

2.24 

.91  1 

.6 

-.0900 

-.176 

.637 

2.28 

.789 

.6  57 

2.33 

.772 

.8 

-.133 

-.250 

.663 

2. 36 

.641 

.661 

2.48 

.594 

1.0 

-.186 

-.333 

.693 

2.50 

.454 

.689 

1.1 

-.216 

-.379 

.710 

2.60 

.334 

.704 

1.2 

-.250 

-.429 

. 728 

2.78 

.170 

. 719 

1.246 

-.267 

-.452 

1 

3.01 

.016 

.726 

^’Reference  (3). 


Tabic  ^ 

COMPARISON  OF  THE  IIIGH-PRANDTI -NUMBER  APPROXIMATION 
TO  NUMERICAL  SOLUTIONS  FOR  IIOWARTH'S  RETARDED  FLOW 
(WATER  WITH  = lO^^F  AND  = 32‘’F; 


Numerlcfll  Solution 

Hlgh- 

Prandtl-Number 

Approximation 

x/1. 

u 

Nu  /Re 

X X 

W 

Nu  /Re^ 

X X 

0.0 

0.923 

0.90f 

.2 

0.890 

0.869 

.4 

1 0.853 

0.820 

.6 

0.803 

0.763 

.8 

0.730 

0.689 

1,0 

0.669 

0.582 

1.1 

0.610 

0.501 

1.2 

0.511 

0.359 

1.246 

0.357 

0.111 

^Reference  (3). 
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V_.  _C0N C LUSIONS 


TIk'  c 1 .iss  i iM  I t sot  hornui  I fonstont -proport  ic-s  Integral  method  of 

Walz  and  Tliwailes  lias  been  extended  to  heated  water  boundary  layers, 

where  t ('mpcr.iture-dependont  viscosity  is  important.  Correlations  of 

the  universal  parameters  K,  li,  and  T as  simple  functions  of  e * and  A 

w 

<illow  computation  of  0,  6*,  and  f"(0)  for  heated  water  boundary  layers 
with  any  surface  shape  and  with  most  temperature  It?vcls  of  interest 
for  underwater  vehicles. 

The  method  has  been  tested  for  the  llowarth  retarded  flow  with 

t^  = 104“F  and  t^  = J2“l'.  With  water,  this  flow  is  ver  / temperature 

dependent;  !u  at  Ing  in  this  flow  eauses  the  separation  point  to  move 

(3) 

It)  percent  further  down  the  plate  ; the  flow  is  highly  nonsir.iilar 
even  without  iR'ating.  The  value  of  0 was  calculated  by  the  integral 
method  to  within  1 percent  nearly  to  the*  separation  point.  The  values 
of  6*  and  f"(0),  both  nearly  exact  at  x = 0,  had  each  attained  an  error 
of  only  2 percent  at  a position  lialfway  to  the  separation  point. 

The  method  shares  several  advantages  and  shortcomings  with  the 
original  metluids  upon  which  it  is  based.  Its  advantages  are  that  the 
momentum  thickness  and  Iieac  transfer  can  be  calculated  accurately, 
easily,  and  economically  for  a wide  variety  of  pressure  gradients,  both 
favorable  and  adverse.  The  shortcomings  are  that  it  may  not  represent 
li  and  T adeciua t e 1 y , particularly  in  regions  of  adverse  gradients,  and 
thus  the  location  of  laminar  separation  may  not  be  determined  accurately. 

The  high-Prandt 1-number  approximation  has  also  been  found  to  be 
accurate  in  computing  the  heat  transfer  in  heated  water  boundary  layers, 
except  in  the  reg  m near  separation. 

Despite  these  limitations,  this  method  extends  a proven  and  re- 
liable approach  to  the  computation  of  laminar  boundary  layers  in  water. 
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Append  lx 

fluid  property  correlations 


The  water  fluid  properties  correlations  used  in  the  numerical  work 
for  this  report  are  those  of  Kaups  and  Smith. The  density  is  taken 
as  constant,  so  only  the  viscosity  ratio,  thermal  conductivity  ratio,  and 
Prandtl  number  correlations  are  needed.  With  T expressed  in  degrees 
Rankine,  and  T^^^  = A91.69°R  (32“F),  these  correlations  are 

y/u^ef  - 1/C33. 15539  - 106.9718715  (T/T^^^)  + 107.7720376  (T/T^^f)^ 

- A0.5953074  + 5.6391948 

= -1.940589  + 5.2220185  (T/T^^^)  - 2.693322  (T/Tj.^c) 

+ 0.4176167  (T/T^gf)^ 

Pr  - 13. 66/[73. 376906  - 208.7474538  (T/T^^^  + 197.7604676  (T/T^g^)^ 

- 68.8626186  (T/T^g^)^  + 7.4779458  (T/Trgf^^^’ 
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